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Cyanopolyynes in hot cores: modelling G305.2+0.2 
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ABSTRACT 

We present results from a time dependent gas phase chemical model of a hot core based 
on the physical conditions of G305.2+0.2. While the cyanopolyyne HC3N has been ob- 
served in hot cores, the longer chained species, HC5N, HC7N, and HC9N have not been 
considered typical hot core species. We present results which show that these species 
can be formed under hot core conditions. We discuss the important chemical reactions 
in this process and, in particular, show that their abundances are linked to the parent 
species acetylene which is evaporated from icy grain mantles. The cyanopolyynes show 
promise as 'chemical clocks' which may aid future observations in determining the age 
of hot core sources. The abundance of the larger cyanopolyynes increase and decrease 
over relatively short time scales, ~10 2 5 years. We present results from a non-LTE sta- 
tistical equilibrium excitation model as a series of density, temperature and column 
density dependent contour plots which show both the line intensities and several line 
ratios. These aid in the interpretation of spectral line data, even when there is limited 
line information available. 
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1 INTRODUCTION 

Hot cores are observed in regions of high mass star forma- 
tion, and are small, hot (~ 200K) and dense (~ 10 6 cm~ 3 ) 



clumps of gas ( e.g.. ISweitzerl|l978l ; iPauls et al. 1 1 1981 for a 
review see, e.g.. Ivan der TakH2003h . The chemistry of these 
regions is unique, characterized by unusually high abun- 
dances of saturated molecules, such as HgO and NH3 (e.j 
iHenkel et al]|l987l : IWalmslev et~allll987l : lOlmi et al.lll993 
This unique chemistry is believed to be a direct result of 
grain mantle evaporation due to an event, such as the switch- 
ing o n of a nearby star or a source within the core itse lf 
(e.g., iKaufman et alj Il998l ; IWvrowski fc Walmslevl Il99rj) . 
During an earlier high density phase of collapse, these 
molecules are formed on the surface of dust grains (e.g., 
IWalmslev fc Schilke|[l99ll ). released into the gas phase dur- 
ing the onset of heating, and are processed into daughter 
species. Since cores are both hot and dense, the time depen- 
dent gas phase chemistry rapidly changes over time scales 
< 100,000 years. Species with time dependent chemistry can 
be used as 'chemical clocks' an d may help in determining 
the a ge of hot core regions (e.g.. iMillar. Macdonald fc Gibbl 
Il997h . 

The origins of some interstellar molecules found in hot 
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cores can be separated into the following three groups. Those 
that are formed in the cold gas: 

e.g., CO, C2H2, C2H4, N a , 2 , CS, 



those formed on grain surfaces: 

e.g., H 2 0, CH4, NH 3 , H 2 S, CH3OH, H 2 CO, C2H5OH, 
CO2, HCOOH, NH2CHO, C2H4, C2H 6 ,H 2 CS, OCS, 



and those formed in the hot gas (e.g.. Ic"harnlev|[l995t ) 

e.g., CH3OCH3, CH2CO, CH3CHO, HCOOC2H5, 
CH3COCH3, (C 2 H 5 )20, CH3OC2H5, C2H5CN, CH3CN, 
HON, SO, SO2, C2H3CN, HC 3 N, CH3NH2, CH 2 NH, SiO. 



Observations of hot cores include molecular line 
and dust continuum emission, and some of the chem- 
ical core models, along with radiative transfer calcu- 
lations, have proven useful in reproducing observed 
column densities (e. g., Millar, Macdonald & Gibb Il997l; 
Kaufman et al. I Il998l : iThompson et al.l 1 19991 : iDotv et al.l 



2002j y Chemical models have been developed that in- 



clude a pre-stellar phase or initial conditions assuming 
the presence of evaporated grain mantle species (e.g., 
iBrown. Charnlev fc Millarll 19881 : iMillar. Herbst fc Charnlevl 
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Millar. Macdonald fc Gibbl 1 19971: IViti fc Williamsl 



1992 



1993 



1999 



Table 1. The initial gas phase abundances, with respect to H nuclei, 
assumed in the chemical model. 
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20041 '). The injection of ice mantle material into the gas 
phase is usually assumed to occur instantaneously since 
it is driven by stellar radiation and the grain tempera- 
ture rapidly increases from 10K (in the cold phase) to 
~ 200K (in the hot core) when the ices evaporate quickly. 
Gas-grain chemical hot core models which assume simple 
molecules like acetylene, C2H2, are released from grain 
mantles, have been used to account for the observed 
abundances of many O and N bearing molecules (e.g., 
ICharnlev. Tielens. fc Millarlll992l ). 

Carbon chain molecules make up a significant fraction 
of the observed interstellar molecules and cyanopolyynes 
HC n N are often detected. The longer chain cyanopolyynes 
HC2n+iN (n=2-5) are associated with regions such as 
dark dust clouds, and are not usually associated with hot 
cores. HCnN ha s been detected in the interstellar medium 
ijBell et al.lll997h . with the shorter chained species, for ex- 
ample HC3N and HC5N, more commonly associated with 
a number of sources including circumstellar shells and star 
forming regions. Cyanodiacetylene, HC5N, is believed to be 
formed inefficiently in hot cores since grain mantles do not 
comprise of larg er unsaturat ed molecul es such as d i acety - 
lene, HCCCCH ljMillarlll997l ). Recently. ISakai et all (|200ct ) 
detected the longer chain cyanopolyynes HC5N, HC7N and 
HC 9 N towards the low-mass pr otostar IRAS 04368+2557 
L1527 . Using a gas grain model, lHassell. Herbst fc Garrodl 
(2008) modelled the chemistry of L1527, and determined 
that the composition of the cloud is most likely a result of 
warm, 30K, chemistry than cold processes. 

The rotational lines from vibrational excited states 
of some cyanopolyynes have been de t ected in hot core 
sources. IWvrowski. Schilke fc Walmsevl (|l999h observed a 
sample of six hot core sources, including G10. 47+003, in 
vibrationally excited lines of HC3N. Observations of the 
Orion Molecular cloud Core (OMC-1) have included the J 
= 12 — » 11 transition of the vq = 1 vibratio nally excited 
state (|Goldsmith. Krotkov fc Snelllll983l . 119851 ). and 5 tran- 
sitions of the vr = 1 vibrationally excited state of HC3N. 
The higher cyanopolyynes are not usually detected, how- 
ever, 10 r otational tran sitions of HC5N and 8 transitions 
of HC7N (|Turnerl Il99l | ) have been observed in the O rion 
core (O MC-1) HC 5 N (|Averv et al.lll979l : lTurnerlll99lh and 
HC7N |Turnerlll99l[ ) have also been detected in the molec- 
ular cloud of Sgr B2. 

Models of the gas phase chemistry are extremely use- 
ful as predictive tools with regards to observations. How- 
ever, extracting molecular abundances from line emission, 
suitable for comparison with the model predictions is not 
a straight forward task. Furthermore, although theoretical 
models can produce abundances for hundreds of molecules, 
observations are limited to specific frequency bands, and 
possibly contain line emission from only a handful of species. 
A broad understanding of the chemistry is thus not of- 
ten possible with limited molecular line data for a specific 
source. 

Modelling tools to aid in the interpretation of spectral 
line data are especially important for large line surveys, and 
efficient methods are needed to obtain physical and chemical 
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parameters, such as column densities, relatively easily even 
with large numbers of species and lines. Currently, collisional 
rate coefficients are available for approximately 24 species 
on the Leiden atom ic and molecular database (LAMDA) 
JSchoier et al.ll2005l). with several others in literature (e.g., 
CH3CN in Greenl 1986h . and so statistical equilibrium cal- 
culations are limited to these. Moreover, observations are 
limited in frequency, and so limit the process further. The 
list of species for which collisional rate coefficients are avail- 
able include CS, OCS, SO, S0 2 , SiO, HCO+, HC 3 N, HCN, 
HNC, H 2 0, OH, CH3OH, and NH 3 . There are many more 
species of astronomical interest that do not appear in the 
LAMDA database, or elsewhere in the literature. 

The aim of this paper is to model the gas phase chem- 
istry of the hot core associated with G305.2+0.2 with partic- 
ular focus on the formation of the cyanopolyynes. Non-LTE 
calculations are used to examine the excitation properties of 
these species. Section[2]outlines the chemical model with the 
cyanopolyyne HC^s^gjN chemistry discussed in Section[3] 
We reveal that the longer chain cyanopolyynes can be pro- 
duced in appreciable abundances assuming the physical con- 
ditions of G305.2+0.2. In Section|3]we present results from a 
statistical equilibrium model which produces integrated line 
intensities. We present a series of density nH2, temperature 
T and column density Nj dependent contour plots which 



will directly aid in the interpretation of observati o ns. N on- 
detections of HC 5 N and HC 7 N from IWalsh et al.l (|2007l ) in 
G305.2+0.2 are explained using the non-LTE model. 



2 CHEMICAL MODEL 

T he physical model for the cor e is similar to that 
of iMillar. Macdonald fc Gibbl (Il997l ) and iThompson et all 
i|l999l ) which assume a spherical core geometry. The den- 
sity and temperature is constant throughout since we want 
to examine the general gas phase chemistries. We as- 
sume a density of nn2 — 8 x 10 4 cm -3 and a temper- 
ature of 200K, which was derived fo r the hot core of 
G305.2+0.2, G305A, in IWalsh et al.1 (|2007T ). with a ra- 
dius of l.lpc at a distance of 3.9kpc. The current mod- 
elling work use s the latest chemical rate coe fficients from 
Woo dall. et al.l i|2007h . The model results o f Walsh ~et al.l 
<|2007h used the earlier iLe Teuff. Millar fc Markwickl (|2000T ) 
version. Also, the current modelling work includes the longer 
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Table 2. Species present in model calculations, grouped into the number n of atoms. 
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chain cyanopolyynes, HC7N and HCgN, along with associ- 
ated species needed for thei r formation. These species were 
not considered previously in IWalsh et al.l (|2007h . 

The ini tial gas phase abundances are given in Ta- 
ble □ (e.g., IMillar, Macdonald fe Gibbl Il997l ; IWalsh et all 
12007]). The initial abund ances are based on previous mod- 
els of hot cores (e.g.. IMillar. Macdonald fc Gibbl 1 19971 . 
iRodgers fc CharnlevlkoO ill Walsh et alj|2007h and the mean 
interstellar values. The source is assumed to fill the beam, 
and so beam dilution does not play a large role in the model 
results. The model contains 245 species, consisting of the 8 
elements H, He, C, N, O, Si, S, and Fe. The reaction rate 
network contains 3184 reactions. Included in the model are 
the cyanopolyynes HC 2n +iN (n=0-4), the carbon chained 
molecules C n H 2 (n=3-9) and the hydrocarbon chains CnH 
(n=l-9). The full set of species are given in Table [2] with 
the initial mantle abundances in Table [1] The species list 



includes many previously observed hot core species along 
with essential species needed for their production. 

Di screpancies between models of dense interstellar 
clouds (jMillar. Leung fc HerbstJI 19871 ') highlight the sensitiv- 
ity of astrochemical models to the reaction rate coefficients 
as well as the reaction network. The reaction rates are taken 
from the most recent UMIST database whic h includes the 
latest laboratory data (|Woodall. et al.l |2007l ). The reaction 
networks of some of the key core species will be discussed 
in more detail below. Many of the rate coefficients are tem- 
perature dependent. There are a select few reactions with 
rate coefficients that are ill-defined outside their prescribed 
temperature range, with rate coefficients exponentially in- 
creasing with decreasing temperature. In these cases the rate 
coefficients are kept constant at the lowest or highest tem- 
perature range value nearest the temperature of interest. 

Numerically, the core is described by a series of 22 con- 
centric shells, with a density and temperature assigned to 
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each. At each depth point or shell, the initial abundances 
are set and the reaction network is calculated over 10 8 years. 
However, the chemistry rapidly evolves over much shorter 
time scales than this. The species abundances are integrated 
along the line of sight and the column density of a species j 
at time t is 



T 1 1 — 



N,(r,e,t) 



rij(s, t)ds 



(1) 



where rij(s,t) is the species abundance and s is the depth 
through the core along the line of sight. The coordinates r 
and 6 specify the position of the pencil beam with respect to 
the centre of the core. We approximate the telescope beam 
using a Gaussian with an area normalized primary beam: 



P(r) = 



27TCT 



2 exp(-^) 



HPBW 



2%/2In2 



(2) 



where r is the radius measured from the centre of the gaus- 
sian and HPBW is the half power beam width. P(r) is nor- 
malized so that 



P(r)rdrd6 = 1. 



(3) 



To evaluate the gaussian beam, column densities are 
firstly evaluated along pencil beams on grid points through 
the chemical core model. These column densities are then 
convolved with the Gaussian, equation ([2]), giving the beam- 
averaged or beam diluted column density Nj(t) 



Nj(t) = 



1 



27TCT 2 



rdrde. (4) 



3 CYANOPOLYYNES 

In a related project, we have begun a survey for HC5N 
in methanol-maser selected hot molecular cores from the 
iPurcell et ai1l2006l list, using the Tidbinbilla 34m telescope 
to measure the 31.95 GHz J=12-ll line. These results will 
be reported separately, when the survey is complete. How- 
ever, we can report that HC5N is frequently observed in the 
cores selected from this list, typically w hen CH3CN 5-4 92 
GHz line emission was also detected bv lPurcell et alj|200d 
Sources have been found to contain both CH3CN and HC5N, 
or neither at all, suggesting their possible use as chemical 
clocks. We will explore the chemical pathways below and we 
show that HC5N may well be a hot core species. We also 
disc uss the likely role of other carbon bearin g molecules. 

iFukuzawa. Qsamura fc Schaefer] (|l998h have demon- 
strated the importance of these neutral-neutral reactions: 



C 2n H 2 + CN -> HC 2n +iN + H, 



(5) 



in the formation of cyanopolyynes in the ISM. Cyanoacety- 
lene, HC3N is often detected in cold gas in molecular clouds, 
and is formed efficiently via the reaction: 

C 2 H 2 + CN^HC 3 N + H, (6) 

with a rate coefficient from ISmith et al.l [2004. Observed 
abundances of C 2 H 2 , (as well as CH3CN and H 2 CO) can- 
not be produced in the hot core gas phase alone, but 
can be accounted for in the cold phase chemistry where 
they are frozen out onto the grains before being re- 
leased into the gas phase (e.g., iBrown. Charnlev fc Millar] 



C2H2 
CN 

HC3N 
HC5N 
HC7N 
HC9N 




Log 10 (t/yrs) 

Figure 1. Time dependent column densities for the chemical 
model of Section \3\ with njj 2 = 8 xl0 4 cm — 3 and a temperature 
of 200K. The different length cyanopolyynes are compared, along 
with the parent species CN and C 2 H 2 . 



ll98StlCharnlev fc Tielenslll992l: Charn lev. Tielens. fc Millar! 
Il992h . lLahuis fc van Dishoeckl (|2000) found high abun- 
dances of hot C 2 H 2 towards massive YSOs with ISO and 
concluded that the most likely origin is the evaporation 
of grain mantles. Maps of C 2 H 2 emis sion at 13.7 microns 
with S pitzer toward Cepheus A East in lSonnentrucker et al.l 
(2007) show that the distribution is likely due to the sputter- 
ing of acetylene in icy grain mantles. Although C 2 H 2 has yet 
to be directly detected in interstellar i ces, the upper limi t 
is high, 10~ 5 with respect to hydrogen dBoulin et al.|[l998l ). 
and the fractional ab undance of C 2 H 2 in Table [T ] is chosen 
to be below this limit. iRodgers fc Charnlev! (2001) have also 
demonstrated with a gas grain model that high abundances 
of C 2 H 2 may be produced in the gas phase provided that 
a large enough abundance of methane is ejected from the 
grain mantles. Inside hot cores, C 2 H 2 is released from the 
grain mantles and reaction (|SJ proceeds relatively easy for 
hot core temperatures (100-300K). 

Results from the time dependent gas phase chemical 
model described in Section [2] are shown in Figure [T] The col- 
umn density of HC3N rises to a peak above N=10 14 ' 5 cm -2 
at t=10 4 ' 3 years, and reaches its steady state abundance af- 
ter 10 6 years. The abundance of HC3N initially increases 
with CN since it is mainly produced via the reaction of 
CN with C 2 H 2 , reaction (JS|. Also plotted are the larger 
cyanopolyynes HCs.r.gN, whose abundances also peak near 
t=10 4,3 years. These species are not usually considered hot 
core species but in Figure Q] they are all reach appreciable 
column densities over N = 10 12 ' 5 cm -2 . The peak column 
densities of the shorter molecules, HC5N and HC7N, are 
one order of magnitude larger than that of HC9N. Also, all 
species reach steady state abundances after 10 6 years. 

In the ISM cyanodiacetylene, HC5N is mainly formed 
in the neutral-neutral reaction 
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CN + H 2 CCCC -* HC 5 N + H. (7) 

The rate coefficient for this re action is 2.55 x 10~ 1 cm 3 
s _1 at a temperature of 300 K (jWoodall. et al.l 120071 ). This 
was confirmed as the dominant reaction in the production of 
HC5N in the hot core results of Figure [T] This reaction has 
not been considered important in the past since it was not 
thought that H2CCCC was formed on the grain surfaces 
or very efficiently in the gas phase. However, H2CCCC is 
formed efficiently via 



C2H + C2H2 



H 2 CCCC + H 



(8) 



with a rate coefficient of 1.5 x 10~ 10 cm 3 s" 1 at 200K. 
The column densities of H2CCCC and C2H are plotted 
in Figure [3] reaching peak column densities of 10 14 ' 8 cm -2 
and 10 14,5 cm -2 , respectively, near 10 4 years. Therefore, the 
abundance of H2CCCC, and thus HC5N, also depend on 
acetylene C2H2. CN is produced in relatively large quanti- 
ties, furthering the efficiency of reaction {7} . CN is produced 
mainly by dissociative recombination of H2NC + , 



H 2 NC + + e" 



CN + H 2 



and at later times after 10 5 years by the reaction 

C + NO -> CN + O. 

H 2 NC+ is formed by 

He + + CO -> O + C+ + He 
C+ + NH 3 -> H 2 NC+ + H, 

and reaction (|10l) is preceded by 

N + OH -> NO + H 
H 3 + +e~ -> OH + H + H. 



(9) 
(10) 

(11) 
(12) 

(13) 
(14) 



Carbon chained polar molecules of the form H2C=C..C 
have been detected in the laboratory and have been identi- 
fied as useful species to radio astronomers be ing stable under 
interstellar conditions. iKillian et al] (|l990h have detected 
propadienlylidene H2CCC and butatrienylidene H2CCCC in 
the laboratory, and suggest that the H2C„ cumulene car- 
benes may be an another important sequence of interstellar 
molecules, comparable with the cyanopolyynes (HC n N ob- 
served in space up to n = 11) and the linear hydrocarbon 
chains C n H (observed up to n=8). H2CCCC is an isomer 
of diacetylene HCCCCH and has bee n detected in the dark 
cloud TMC-1 jKawaguchi et all 1991 ) and th e circumstellar 
shell of IRC +10216 dCernicharo et aflll99ll ). 

The increase of the abundance of HC7N in Figure [T] is 
due to the reaction network 



C2H2 + C4H 
CN + C 6 H 2 



C 6 H 2 + H 
HC 7 N + H. 



(15) 
(16) 



Once again the abundance of the cyanopolyyne HC7N is 
directly related to acetylene. Similarly, the increase in the 
abundance of HC9N can be accounted for by the reactions 



CsH2 + C2H 
CN + C 8 H 2 



C8H2 + H 
HC 9 N + H, 



(17) 
(18) 



and this can be related to the HC7N reaction network 
(reaction I15p via the C6H2 molecule. The column densi- 
ties of the molecules responsible for the production of the 
cyanopolyynes are plotted in Figure [3] and are compared 
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Figure 2. As for Figure [l] but for a comparison of the higher 
order cyanopolyynes HC n N (n>3) with their parent species. 



with the daughter species HC5N, HC7N and HC9N in Fig- 
ure The abundance of H2CCCC increases to a peak of 
10 14 ' 8 cm -2 at 10 4 years. The abundance of the daughter 
species HC5N also increases reaching a peak of 10 14,2 cm~ 2 
near 10 5 years. The peak in HC5N occurs at a later time then 
its parent H2CCCC. Similar relationships between CeH2 and 
HC7N, and between CsH2 and HCgN, can also be seen. 

The formation of the cyanopolyynes, HC5N , HC7N and 
HC9N, are sensitive to the abundances of CN and the car- 
bon chain molecules H2CCCC, C6H2, and CsH2. Further- 
more, the abundance of the carbon chain molecules are de- 
pendent on acetylene C2H2. Thus one would expect, that 
if C2H2 is present in the gas then the larger cyanopolyynes 
may also be present. Figure [4] illustrates this, showing the 
resulting abundances when no C2H2 is initially evaporated 
from grain mantles. The acetylene abundance initially in- 
creases to a column density ~ 10 14 cm -2 near 10 4 ' 3 years. 
HC 3 N, HC 5 N, HC7N and HC 9 N reach peak densities of 
10 135 cm- 2 , 10 12 ' 7 cm- 2 , 10 118 cm- 2 , and 10 104 c m - 2 , re- 
spectively. These column densities are a magnitude (or in 
some cases less) smaller then those of Figure [T] 

Since CH3CN emission is often used as a signature of 
the presence of a hot molecular core, we have investigated 
its relative abundance compared to HC5N as a function of 
time. CH3CN is efficiently produced via 



CH 3 CNH + + e" 



CH3CN + H 



(19) 



and the abundance of CH3CNH 
sociation: 

CRt + HCN -> CH 3 CNH + + hv. 



is driven by radiative as- 



(20) 



The methyl cyanide chemistry is shown in Figure [5] For 



the first 10 years, CH3 is produced by 
H 3 + + CH3OH -> CUj + H 2 + H 2 



(21) 
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Figure 3. As for Figure[T] but for the parent species of the higher 
order cyanopolyynes HC n N (n>3). 
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Figure 5. As for Figure [T] but shows the parent species respon- 
sible for the formation of CH3CN. 




after 10 5 years. These are preceded by 



Log 10 (t/yrs) 

Figure 4. As for Figure[T] however, C2H2 is assumed to form in 
the gas only, with no initial abundance related to grain mantle 
evaporation. 



H 2 + HCN+ 
C+ + NH 3 
C+ + NH 2 



HCNH+ + H 
HCN+ + H 2 
HCN+ + H. 



(25) 
(26) 
(27) 



In Figure [5] the destruction of the mantle species NH3 
and CH3OH begins near 10 4 years, signalling the start of the 
rapid gas phase chemistry between 10 4 and 10 5 years. As re- 



action (|23[) proceeds, the HCN abundance increases. CH3 is 
strongly dependent on CH3OH which is highly abundant 
before 10 4 years and so remains in the gas phase with a 
constant column density of 10 12 cm -2 between 10 2 and 10 4 
years. When the abundance of HCN increases, reaction (|20|l 
proceeds and the abundance of CH3CNH 4 " increases, allow- 
ing the formation of CH3CN via reaction (|19[) . 

The abundances of HC5N and CH3CN are not directly 
linked via their respective reaction networks, however the 
abundances of both are strongly dependent on their parent 
grain mantle species. The abundances of both increase at 
the onset of the rapid gas phase chemistry caused by the 
destruction of the grain mantle species. As a useful com- 
parison, we plot the time dependent HC5N/CH3CN ratio in 
Figure [S] The ratio increases considerably between 10 4 and 
10 years and so could prove useful in determining the age 
of hot cores within this range. 



and 

H 2 + CH 2 f -» CH^+H (22) 
thereafter. The HCN abundance is driven by 
NH 3 + HCNH + -> HCN + NHj (23) 
for the first 10 years and 

HCNH + + e" -f HCN + H (24) 



4 MOLECULAR LINE EXCITATION MODEL 

In this section, a molecular excitation model is used to cal- 
culate the excitation of HC 3 N, HC 5 N and HC 7 N. The col- 
lisional rates (per second pe r molecule) with H 2 are taken 
from the LAMDA database (|Schoier et al.ll2005l ). Only col- 
lision rates for the ground vibrational state are currently 
available in the database and it outside the scope of this 
work to consider the higher vibrational transitions. We note 
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Figure 7. Line intensities and ratios of HC 3 N for the J=9-8 (81.88GHz), J=10-9 (90.98GHz) and J=ll-10 (100.08GHz) lines calculated 
by the excitation model. The first row shows the line intensities in units of K km/s as a function of temperature and column density 
(normalised by line width) assuming a density of njj2 = 8 X 10 4 cm — 3 . In the second row, three line ratios as a function of temperature 
and column density for constant density 8 xl0 4 cm — 3 are plotted. In the third and fourth rows, the corresponding plots as a function of 
density for a constant temperature of 200K are shown. 



however, that although highe r vibrational lines of HC3N 
have been observed in Sg B2 (iGoldsmith. Krotkov fc Snell 
1 19831 . 1 19851 ). I Walsh et all (|2007h did not detect any higher 
vibrational transitions of HC 3 N, HC 5 N or HC 7 N in G305A. 
Here we only consider the ground vibrational states and will 



have to wait until the collision rates are available for higher 
states to update our model in the future. 

The collision rates depend on the cross section of HC3N 
and the density of H2. The LAMDA database also contains 
the relevant transitional frequencies, rotational levels and 
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Figure 8. As for FigureEJbut for the J=9-8 (23.96GHz), J=10-9 (26.63GHz) and J=ll-10 (29.29GHz) lines of HC 5 N. 



Einstein A-coefiicients for 24 molecules compiled from 
sources such as the JPL Molecular Spectrosc opy database 
JPickett et al.lll998l ) and the CDMS catalogue (|Muller et all 
I2005I ). The escape probability method is used here to sim- 
plify the statistical equilibrium equations (see for example 
ISchoier et ai1l2005h . with the model outputting the inten- 
sity of a line for a given density and temperature. Similar 
model calculations for ot h er spe ci es of interest can be found 
in, for example, iLisztl l|2006l ). iLique. Cernicharo &; 



(2006|V iDaniel. Cernicharno fc Dubernetl (1 20061). 

Ivan der Tak fc Hogerheiidel (|2007h . and IWalsh et all 
( 200 it ). 



In Figure [7] line intensities and line ratios are plotted 
for the 81.88, 90.98 and 1 00.08GHz rotationa l lines of HC 3 N. 
The LAMDA database (|Schoier et al.ll2005l ) contains colli- 
sion rates of the first 21 levels (9-181GHz) of HC 3 N (how- 
ever the higher transitions are affected by truncation) . These 
plots are useful in interpreting observations since a measured 
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Figure 6. Time dependent HC5N/CH3CN ratio for model in 
Figure 1. 



line intensity, or line ratio can directly be compared on the 
plot, and the column density and/or density may be de- 
termined. In the first row, the individual line intensities in 
units K kms -1 are given for a density of nn2 = 8 x 10 4 cm -3 
as a function of temperature and column density. Each plot 
shows the same overall trends, with the intensities becoming 
relatively constant with temperature for given normalised 
column densities below 10 14 cm -2 . A single observed line in- 
tensity could be directly compared to these, after telescope 
beam dilution effects h ave been taken into account. As illus- 
trated in, for example, (I Walsh et al.ll2007l ). plotting multiple 
lines from non-LTE models allows for the determination of 
the column density, as well as putting limits on the density 
and temperature of the region. 

We also show in the second row of Figure [7] the line 
ratios as a function of temperature for nH2 = 8 xl0 4 cm~ 3 , 
which do not require any correction for beam dilution ef- 
fects. The intensity of the 100.08GHz (J=ll-10) line is 
strongest (for a given temperature and density), then the 
90.98GHz (J=10-9) and 81.88GHz (J=9-8) lines. In the hot 
core regime, for temperatures of the order 150K and above, 
the ratios are relatively constant with column density and 
so potentially allowing one to read off the column density 
directly using an observed line ratio value. 

In the third row of Figure [7] the line intensities are 
shown again, but this time as a function of density for 
T=200K. Each contour becomes constant with column den- 
sity as the density increases, and the LTE approximation 
is approached with increasing nu2- One could safely assume 
the LTE approximation, and use a rotation diagram, for den- 
sities above 10 6 cm -3 and if the line width normalised HC3N 
column density is expected to be in the range 10 12 cm -2 to 
10 16 cm~ 2 . Finally, in the bottom row of Figure [7] the line 
ratios are plotted as a function of density. As the column 
density decreases, the contours become constant with den- 
sity. The trend towards LTE with increasing density is also 
clear. 

Depending on what information is available, either se- 



ries (row) of plots in Figure [Jj may be the most useful in 
interpreting observations. If either the temperature or den- 
sity is known than the complexity is reduced greatly, and 
the corresponding plots with the appropriate temperature or 
density can be produced. Notwithstanding, however, even if 
these are unknown, the combination of all four plots can 
still be useful in determining limits on the density, tem- 
perature and column density. Collision rates for the larger 
cyanopolyynes HC5N, HC7N, HC9N etc. have yet to be pub- 
lished. To calculate the exc itation of HC5N, we have used the 
collisional rates of HC3N jSchoier et al.l [20051 ). however the 
cross section for HC5N is scaled by a fact or of 1.5 to account 
for the greater geometric length (e.g., lAverv et al.l Il979l ; 
ISnell et a.1.1 Il98lh . The energy levels and Einstein A coeffi- 
cients for HC5N ha ve been taken from the CDMS database 
|Muller et alj|2005h . The collision rates for HC 3 N are avail- 
able for the first 21 lines and in the case of HC5N this in- 
cludes all transitions from 2GHz to 54GHz. In Figure [S] the 
corresponding line intensities and line ratios are plotted for 
the 23.96, 26.63 and 29.29GHz lines of HC 5 N, showing sim- 
ilar characteristics as for HC3N. A similar calculation has 
also been performed for HC7N, after scaling the collision 
rates of HC3N by 2 to account for the larger cross section 
of HC7N (plots not shown here). These plots show similar 

features as those for HC3N and H C5N. 

In the ATCA observations of IWalsh et""ai1 <|2007f ) the 
HC 5 N(7-6) 18.64GHz and HC 7 N(15-14) 16.92GHz lines 
were not detected in G305A, although Figure[JJsuggests that 
these species are likely produced in this region. We will now 
address this issue. The ler upper limits of the integrated 
intensity of these non-detections from the ATCA observa- 
tions were <0.11Kkm/s and <0.03Kkm/s, for HC 5 N(7-6) 
and HCtN(15-14), respectively. Using the non-LTE model, 
we plot the upper limit intensities of each of these lines 
in Figure [9] for a temperature of 200K, and varying den- 
sity. We also plot the observed HC 3 N(2-1) 18.20GHz line 
for comparison, this line was al r eady considered in the 
non-LTE analysis of IWalsh et al] (2007). For a density of 
nH2 = 8 x 10 4 cm -3 , the following limits are imposed on the 
normalized column densitie s; NhCrn/^U < 3.1 6 x 10 12 and 
N H c 7 n/5V < 3.89 x 10 11 . In lWalsh et all (|2007h an age range 
for G305A of 2 x 10 4 < t core < 1.5 x 10 5 year was imposed, 
based on all the detected lines. From Figure[TJ this age range 
implies column densities of 10 13 cm -2 < Nhc 5 n < 10 14 cm -2 
and 2.5 x 10 12 cm -2 < N H c 7 n < 6.3 x 10 13 cm~ 2 , with the 
lower limits defined by the later core age of 1.5 x 10 5 years. 
After, assuming a line width of > 3.2km/s for the lines 
in Figure [5] which is not unreasonable, the column den- 
sities of HC5N and HC7N agree with Figure [1] for a core 
ra nge of the order 1. 5 x 10 5 years. The NH3 observations 
of IWalsh et all i|2007f ) favor a later age for G305A of the 
order 10 5 years, also. Around this time the cyanopolyynes 
are being destroyed, and their column densities decrease 
rapidly. This could explain why these species were not de- 
tected in lWalsh et al.l (|2007T l. since there abundance was not 
high enough to be detected. 



5 DISCUSSION 

The results from the chemical model in Sections [2] and [3] 
show that the cyanopolyynes switch on after 10 2 years and 
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Figure 9. Comparison of cyanopolyynes from (Walsh et al .120071) 
using a non-LTE excitation model. 



peak around 10 4 ' 5 years, before being rapidly destroyed be- 
fore 10 5 ' 5 years. Acetylene was shown to be a key species 
in the formation of the cumulene carbenes HbCn, which re- 
act with CN to form the cyanopolyynes. Since acetylene is 
abundant in large quantities, as it has been evaporated from 
grain mantles, and CN is formed easily in the gas phase, 
the cyanopolyynes should also be found in hot cores. The 
formation and subsequent destruction of the longer chain 
cyanopolyynes are relatively short lived, which may explain 
why these are not often observed in hot core sources. How- 
ever, if they are observed, they may prove very useful as 
chemical clocks, as they could help pinpoint the age of the 
hot core. Their detection in sources with similar physical 
conditions to those of G305A in Section [2] should imply ages 
of the order of 10 4 ' 5 years. 

The calculated time dependent chemistry in Section 
[2] are for specific physical conditions; iih2 = 8 x 10 4 cm -3 
and T=200K. We now examine the sensitivity of the chem- 
ical model results to changes in density and temperature so 
that these results may be extended to hot cores with dif- 
ferent physical conditions to G305A. Since the conditions 
inside the core are constant throughout, and the source fills 
the beam and beam dilution effects are not so important, 
changes in the size of the core, will only scale the column 
densities and not change the time dependent features. Since 
the reaction rates are temperature dependent, a change in 
temperature will effect the model results. The reaction rates 
for two body reactions are of the form 



k = a(T/300)' 3 exp(7/T) 



(28) 



where a, B and 7 are constan ts taken from the UMIST 
database (|Woodall. et al.ll2007l ). The reaction rates for re- 
action ([5]), are the same regardless of the length of each 
cyanopolyyne but vary with temperature. A reduction in 
temperature to 100K, increases the reaction rate by 1.6 and 
therefore increases the abundances. This is shown in Figure 
[ES1 

Since we have only considered the physical conditions 
for G305.2+0.2, we now increase the density to illustrate 
that longer chain cyanopolyynes should also be formed in 
higher density cores. In Figure [TT] we consider an increase 
the density by a factor of 10 to nH2 = 8 x 10 5 cm -3 . The 
species abundances are enhanced by at least a factor of 10. 
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Figure 10. Comparing the time dependent chemistry for model 
with n H 2 = 8 xl0 4 cm~ 3 and temperatures of 100K and 200K. 
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Figure 11. Comparing the time dependent chemistry for model 
with a temperature of 200K and densities of njj2 = 8 xl0 4 cm -3 
and n H 2 = 8 xl0 5 cm -3 . 



The time dependent features follow the same movements as 
for the lower density G305A case. 



6 SUMMARY 

In the core chemical model presented here, with density 
nH2 = 8 x 10 4 cm -3 and a temperature of 200K, the 
cyanopolyynes show highly time dependent chemistry be- 
tween 10 2 and 10 6 years. Although the larger cyanopolyynes 



Cyanopolyynes in hot cores: modelling G305.2+0.2 11 



HC n N (n> 3) are usually not observed and have not been 
considered in theoretical models in the past, we have shown 
that they can be produced under hot core conditions. Their 
abundance is directly linked to acetylene, which is highly 
abundant in the earlier stages of the core evolution when 
it is evaporated from grain mantles. The reactions respon- 
sible for processing acetylene into the daughter species 
cyanopolyynes are efficient under the physical conditions 
of G305A, as well as the higher density and temperature 
cases considered in Section [5] The abundance of the larger 
cyanopolyynes increase and decrease over a short period of 
the order of 10 2 ' 5 years, and are relatively short lived com- 
pared with other more commonly observed species. This may 
explain why they have been observed in so few sources. Since 
they are only present in the gas phase for a relatively short 
time, they make useful species as chemical clocks, if the 
source age lies within this window. HC3N, however, is pro- 
duced in large abundances much earlier on than the larger 
cyanopolyynes which may explain why it has been detected 
in several sources. 

Through the use of a molecular excitation model, we 
presented a series of line intensity and line ratio plots, use- 
ful for the interpretation of spectral line data of HC3N and 
HC5N. Given a series of line intensities measured from a 
source, these and the associated line ratios can be used to 
look up the column density of a species. Although there wer e 
non-detections of HC 5 N and HC 7 N in IWalsh etafl (|2007h . 
we show that this is most probably due to these species 
having low column densities which are below the observing 
range of the ACTA. 
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